Massless Dirac-Weyl Fermions in a T 3 Optical Lattice 
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We propose an experimental setup for the observation of quasi-relativist ic massless Fermions. It 
is based on a T 3 optical lattice, realized by three pairs of counter-propagating lasers, filled with 
fermionic cold atoms. We show that in the long wavelength approximation the 75 Hamiltonian 
generalizes the Dirac-Weyl Hamiltonian for the honeycomb lattice, however, with a larger value of 
the pseudo-spin S = 1. In addition to the Dirac cones, the spectrum includes a dispersionless branch 
of localized states producing a finite jump in the atomic density. Furthermore, implications for the 
Landau levels are discussed. 
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In the past decade, ultra cold atoms have emerged as 
a fascinating new area linking quantum optics with solid 
state physics [1]. Essentially, these are the only quantum 
many-body systems for which the particle interaction is 
both rather precisely known and controllable. In particu- 
lar, cold atoms confined in optical lattices (OLs) [2] often 
present systems with crystalline structure in various spa- 
tial dimensions d = 1,2,3 described by textbook models 
from solid state physics with tunable parameters. This 
implements Feynman's pioneering idea of quantum sim- 
ulations using one physical system to investigate another 
one [3]. A celebrated example [4] is the optical realiza- 
tion of the Mott transition, a well-known phenomenon in 
solid state physics, describing the transition from a metal 
to an insulator with increasing interaction strength. Fur- 
thermore, the possibility to realize an effective magnetic 
field by rotation of cold atoms in OLs [5] has opened up 
prospects of studying other fundamental phenomena in 
a controlled manner such as the fractional quantum Hall 
effect in d = 2 [6]. 

The recent preparation of single layers of graphene [7] 
has attracted considerable attention, since this solid state 
system displays quasi-relativistic motion of electrons on 
a two-dimensional honeycomb lattice (HCL). However, 
e.g. due to disorder or impurities, many properties of 
real graphene cannot fully be accounted for by the ideal- 
ized Dirac-Weyl Hamiltonian. In this Letter we present 
a detailed study of the T3 lattice [8] and show that 
cold fermionic atoms in such an OL indeed behave as 
quasi-relativistic massless Dirac-Weyl Fermions. Yet, the 
% lattice replaces the pseudo-spin S = 1/2 of Dirac-Weyl 
particles in the HCL by the larger value 5 = 1. As one 
of its crucial features, the T3 lattice exhibits nodes with 
unequal connectivity. The corresponding class of two- 
dimensional lattices, specifically bipartite lattices, has 
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been studied extensively in the past, with a particular fo- 
cus on topological localization [8, 9], frustration in a mag- 
netic field [10, 11], and effects of spin-orbit coupling [12]. 
The T3 lattice, illustrated in Fig. la, has a unit cell with 
three different lattice sites, one six- fold coordinated site 
H, called hub, and two three-fold coordinated sites A and 
B, called rims. All nearest-neighbor pairs are formed by 
a rim and a hub. The energy spectrum [8, 10] of the 
T% lattice exhibits particle-hole symmetry and is char- 
acterized by three branches. Two of them are linearly 
dispersing near the if -points, in direct analogy to the 
energy dispersion known from the HCL [7]. The third 
branch is dispersionless, constantly equal to zero energy. 
These localized states correspond to non-trivial solutions 
of the Bloch equations for which the wave functions van- 
ish identically on the hub sites and have opposite ampli- 
tudes on the two different rim sites. Localization giving 
rise to these dispersionless states has a purely topological 
origin [9] and is quite robust against disorder [8]. 

In this Letter we focus on low energy, long wavelength 
properties of non-interacting Fermions in the T3 lattice 
near half filling. As a main result we find that in this 
parameter range the dynamics is governed by a Hamilto- 
nian of the Dirac-Weyl form 

H = v F S- (p- ^A) . (1) 

Here, v F is the Fermi velocity, p = —i(d x ,d y ,Q) the mo- 
mentum operator (we set h = 1) in the lattice £?/-plane, 
and A the vector potential associated with a magnetic 
field B = (d x A y — d y A x )e z perpendicular to the plane. 
Equation (1) strikingly resembles the Hamiltonian for 
Fermions in graphene. However, as decisive difference, 
the pseudo-spin vector S = (S x ,S y ,S z ) now has total 
spin S = 1 reflecting the three inequivalent lattice sites 
per unit cell in T3 , compared to two in the HCL. In the 
basis of the eigenstates of S z the spin operators of (1) 
can be expressed as 3 x 3 matrices 



2 




:o:c 



E + (k) 



FIG. 1: (Color online), (a) The %, lattice. It is char- 
acterized by translation vectors vi = (3/2; — y/3/2)£o and 
V2 = (3/2; v / 3/2)^o, with lattice constant £$. Open circles 
mark the two sublattices A and B, forming a HCL. Solid cir- 
cles mark the hub sites H forming a (larger) triangular lattice, 
(b) Contour plot of i£+(k), cf. Eq. (3b). The dashed hexagon 
defines the first Brillouin zone, K = 2tt£ " 1 (1/3; -y/3/9) and 
K' = 27t£ " 1 (1/3; V3/9) are two non- equivalent Dirac points. 
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which satisfy angular momentum commutation relations. 

For zero magnetic field the Hamiltonian (1) is obtained 
by starting from the Schrodinger equation for the T3 lat- 
tice in the tight-binding approximation: 

Mh(Rh) = -t^^ A (RH+T J ) + ^ B (RH+r J +i), (2a) 
3 

M a (R a ) = -t^f H (Ra-r j ) J OL G {A,B}. (2b) 

j 

Here, $ Q (R a ) is the amplitude of the wave function on 
sublattice a = A,H,B, and the Tj connect nearest neigh- 
bors. Solving Eq. (2), we obtain the energy spectrum 



E (k) =0 

E±(k) =±t[6 + 4{cos[(v 2 
+ cos[vi • k] - 



Vl)-k] 

cos[v 2 -k]}] 1/2 



(3a) 



(3b) 



where E± exhibit a linear dispersion, with Fermi velocity 
v~p = 3t£o/y/2 [13], about either of the Dirac-points K 
and K', as seen in Fig. 1(b). Next we perform a long 
wavelength approximation about K 



* a (R a )=e lK ' K «i(R a ) 



(4) 



by separating out the slowly varying part i(; a of the am- 
plitude. Similarly, we can treat the neighborhood of K'. 
In the absence of short wavelength scattering processes, 
wave functions from the vicinities of K and K' can be 
mapped onto each other by interchanging the two rim 
components so that we focus on the vicinity of K in the 



following. Inserting ty a into Eq. (2) finally yields the re- 
sult (1) when arranging the components a = {A,H,B} 
into a pseudo-spin triplet and expanding to linear order 
in k = |k - K| for k <C |K|. 

Since [H, J z \ = 0, where J z = L z + S z is 
the z-component of the total angular momentum and 



-idp, we can choose simultaneous eigenstates of 



H and J z . At non-zero energy E 
polar coordinates (r, ip) 



zL^f^ these read in 
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with integer m, the eigenvalue of J z . Here, J m (x) is a 
Bessel function, £± = sgn(m ± + ), where the infinitesi- 
mal + ensures the correct values for m = 0, and N is a 
normalization constant. The eigenstates at zero energy 
are given by 






v e+J| m |+^(^)e i (-+ 1 )^ 



(6) 



where M is again a normalization constant. This solution 
has a finite amplitude only on the rim sites which thus 
are topologically disconnected. The zero energy solution 
(6) is infinitely degenerate with respect to the quantum 
numbers n and m. The properties of this level, e.g. the 
delta-like singularity in the density of states, are distinc- 
tive features of the T% lattice and its Dirac- Weyl behavior, 
as we shall see below. 

Let us now consider the Hamiltonian (1) in a perpen- 
dicular, effective magnetic field. We use the symmetric 
gauge A = B{— y, 0)/2. Denning the magnetic length 
£b = y /f 2c/eB, the cyclotron frequency uo c = vy/^b, 
and the dimensionless radial coordinate p = r/£#, the 
Schrodinger equation for eigenenergy E = uo c £ reads 




0. 



(7) 



Here A = ^e~ 1(p (~d p + p - L z /p). The set of first 
order differential Eqs. (7) can be recast into second- 
order differential equations. For the hub component 

^h(/w) <x 0H(p)e im< ^ we get 
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-T+^ 2 UhW=0 (8) 



and similar equations hold for 4>a/b(p)- The general so- 
lution at e 7^ is given by 




/A/a pW-t-L^-ip 2 ^™- 1 ^ 

A/" H sgn(£) p\ m \l)™\p 2 )e™v 
\Af B pl m l + «+4™'^ + (p 2 )e i(m+1) 



(9) 
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with integer m and positive integer n. The J\T a are appro- 
priate numerical factors [14]. Here, L b a (x) is an associated 
Laguerre polynomial of order a > 0, and n± = #(m±0 + ), 
where is the unit-step function. The eigenenergies of 
the states (9), the Landau levels, are found to read 



n + \m\0(— m) + - . 



(10) 



They scale with the square root of the index n and m 
and display a finite energy value even for n = m = 0. 
These peculiarities are in analogy with the unusual Lan- 
dau levels of the HCL that have received broad attention 
in the context of graphene [15, 16]. 

Equation (7) deserves special care at zero energy. 
Then, two solutions are allowed: a zero-energy Landau 
level and a topological solution. The former is character- 
ized by a single non- vanishing component 
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FIG. 2: (Color online) (a) Distribution of the laser field inten- 
sity for generating a T3 lattice, (b) Cut of the field intensity 
along the y-axis. 
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with m > 1 [17]. The other solution has a vanishing hub 
component and reads 





\Sr B p |m|+ ^4"'^ + (p 2 )e i(m+1 



e 2 



(12) 



with n > 1 [18]. The zero-energy Landau level (11) has 
the same degeneracy as the other Landau levels at finite 
energies (10). In the case of the HCL a similar solu- 
tion is responsible for the atypical quantum-Hall effect 
observed in transport experiments through graphene lay- 
ers [19]. The Atiyah-Singer index theorem [20] connects 
these zero-energy Landau levels with the structure of the 
Dirac-Weyl Hamiltonian. Thus, contrary to the finite en- 
ergy Landau levels, they are robust against fluctuations 
of the electrostatic and magnetic fields. On the other 
hand, the existence of the topological solution (12) is 
due to the bipartite structure of the % lattice. The two 
types of zero-energy solutions cannot be distinguished on 
the mesoscopic scale. Their level degeneracies just add 
up and enhance the total degeneracy at zero energy. This 
feature should be observable in quantum Hall effect ex- 
periments. 

Experimentally, the T3 structure can be realized as an 
OL created by three counter propagating pairs of laser 
beams with the same wavelength A = 3/2^o which divide 
the plane into six sectors of 60°. All six laser beams are 
linearly polarized with the electrical field in the xy-pl&ne. 
Orienting the polarization of one pair of lasers along the 
y axis, Ei = (0,^,0), the other two pairs are obtained 
by rotating E by 120° around the z-axis. This arrange- 
ment produces the profile of laser intensity depicted in 
Fig. 2(a); Fig. 2(b) shows a cut of the field intensity 
along the y-axis. It is apparent that the hopping proba- 
bility between rim sites is exponentially small compared 
to the hopping probability between neighboring hub-rim 



pairs. Contrary to other theoretical proposals for the ob- 
servation of massless Fermions in OLs, we do not require 
the presence of time-dependent potentials [21], staggered 
gauge fields [22], external gauge fields [23], and the pro- 
posed laser configuration does not create spurious lattice 
minima [24]. 

The tight-binding Hamiltonian and its long wave 
length approximation (1) are valid for a T3 lattice popu- 
lated by single-component fermionic atoms, e.g. 40 K or 
6 Li. Indeed, for single-component Fermions, the atomic 
collisions are negligible at low temperature [1]. From 
the experimental point of view, time-of-flight imaging via 
light absorption [25] can be used in order to detect the 
presence of massless Dirac-Weyl Fermions. The harmonic 
trap potential V(r) = mco 2 r 2 /2 confining the fermionic 
cold atom gas is ramped down slowly enough for the 
atoms to stay adiabatically in the lowest band while their 
quasi-momentum is approximatively conserved. Under 
these conditions, free Fermions expand with ballistic mo- 
tion and, from the measured absorption images, it is pos- 
sible [26, 27] to reconstruct the initial reciprocal-space 
density profile of the trapped gas. Then, the local den- 
sity approximation is typically well satisfied and the local 
chemical potential can be assumed to vary with the radial 
coordinate as /i(r) = /io — V(r), where /iq is the chemical 
potential at the center of the trap. For a system of cold 
atoms at temperature T, the atomic density is uniquely 
determined by the chemical potential 



n(/i) 



fjb) dk . 



(13) 



Here So is the area of the first Brillouin zone of the % lat- 
tice, and /(k,/i) = [exp[(£ a (k) - fi)/k B T] + l]" 1 is the 
Fermi distribution function, where E a (k) is the energy 
spectrum of the T3 lattice, cf. Eq. (3). Figure 3 shows 
the atomic density n as a function of the chemical po- 
tential \i. The contribution from the highly degenerate 
topological band (6) manifests itself at \i = as a sharp 
jump in the atomic density. This feature is specific to 




FIG. 3: (Color online) The number of atoms n per unit cell at 
zero temperature as a function of ji/t, ignoring physical spin 
and valley degeneracies, where t is the nearest neighbor hop- 
ping amplitude, cf. Eq. (2). Inset: density of states, dn/dp, 
versus fi/t. 

the % lattice. For small finite /i we see that n increases 
(decreases) proportional to /i 2 , which reflects the linear 
dispersion of massless Fermions near the band center as 
well as particle- hole symmetry [28] . On the contrary, for 
values of \i close to the maximum or minimum of the en- 
ergy band, i.e. far away from the band center, where the 
long wave length approximation can no longer be applied, 



n varies proportional to ji. 

In conclusion, we have discussed an experimental setup 
for the observation of Dirac-Weyl Fermions in a 73 OL. 
In particular, we have shown that in the low energy and 
long wavelength approximation the %, Hamiltonian de- 
scribes massless Dirac-Weyl Fermions. This generalizes 
results known for the case of the HCL, however, with 
pseudo-spin S = 1. The Dirac cones manifest themselves 
in a quadratic dependance of the atomic density on the 
chemical potential. Besides the Dirac-Weyl Fermions, 
the Ts energy spectrum features a dispersionless branch 
of localized states. These states give rise to a pronounced 
finite jump in the atomic density, a hallmark of the T3 lat- 
tice. We have also studied the effects of a effective mag- 
netic field perpendicular to the lattice. This leads to 
the formation of unusual Landau levels resembling those 
known for the case of the HCL. In particular a zero- 
energy Landau level coexists with the dispersionless level 
of localized states. 
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